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摘 要 :为 了 精确 预测 民用 机 场 跑道 的 剩余 使 用 寿命 ,利用 数据 融合 技术 将 两 种 数据 集 进 行 联合 分 
析 , 采 用 贝 叶 斯 概率 预测 方法 对 机 场 跑道 定期 检测 数据 进行 更 新 。 考 虑 机 场 道 面 退化 过 程 的 随机 
性 和 动态 性 ,建立 了 一 种 机 场 道 面 性 能 退化 的 动态 半 马 尔 可 夫 随 机 过 程 模型 ,利用 生存 分 析 对 半 马 
尔 可 夫 过 程 模型 的 转移 概率 进行 了 估计 。 考 虑 飞行 交通 量 和 道 面 厚度 的 双重 影响 ,采用 某 地 方 民 
用 机 场 2007 年 ~2017 年 的 道 面 性 能 定期 检测 数据 ,分 析 了 两 种 影响 因素 的 作用 。 利 用 半 马 尔 可 夫 
过 程 模型 的 转移 概率 对 某 地 方 民 用 机 场 跑道 性 能 退化 过 程 进行 了 预测 ,并 采用 马尔 可 夫 链 蒙特 卡 
洛 ( Markov Chain Monte Carlo, MCMC ) 模拟 技术 ,基于 贝 叶 斯 统计 分 析 方 法 ,利用 不 定期 检测 数据 得 
到 的 先 验 信 息 对 半 马 尔 可 夫 过 程 模 型 的 转移 概率 进行 了 更 新 ,将 更 新 后 的 模型 应 用 于 民用 机 场 跑 
道 性 能 预测 ,并 将 预测 结果 和 未 更 新 的 动态 半 马 尔 可 夫 过 程 模型 预测 结果 进行 了 对 比分 析 。 结 果 
表明 ,基于 MCMC 的 贝 叶 斯 分 析 方法 能 够 融合 定期 和 不 定期 检测 数据 确定 动态 半 马 尔 可 夫 过 程 模 
型 的 参数 先 验 信 息 , 可 以 有 效 地 更 新 机 场 跑道 的 性 能 预 估 模型 ,提高 模型 的 预测 精度 。 
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Abstract :In order to accurately predict residual service life of civil airport runway ,two data sets were ana- 
lyzed by using the data fusion technology and joint estimation approach. The periodic test data of airport 


runways were Updated with Bayesian probability prediction method with the irregularly nondestructive 


收 稿 日 期 :2020-10-02 修 回 日 期 :2021-12-06 

基金 项 目 :河南 省 高 等 学 校 重 点 科研 项 目 ( No.21B580008 ) ;河南 省 科技 攻关 项 目 (No. 182102310747 ) ;郑州 新 郑 国际 机 场 三 期 扩建 工程 技术 
研究 项 目 (No. NFXQGZ-033) 

通信 作者 : 郭 成 超 ,教授 。E-mail: guocc@ zzu. edu. cn 

引用 格式 : 魏 保 立 , 郭 成 超 , 王 复明 ,等 . 机 场 跑道 性 能 退化 随机 过 程 模型 的 贝 叶 斯 数据 更 新 方法 [可 . 应 用 力学 学 报 ,2022 ,39(2) :356-366. 
WEI Baoli ,CUO Chengchao, WANG Fuming,et al. Data updating method with Bayesian analysis for stochastic process model of performance 
degradation of airport runway[ J |]. Chinese journal of applied mechanics ,2022 ,39(2 ) :356-366. 


testing data of runways. Considering the randomness and dynamics of airport pavement degradation, a 
dynamic semi-Markov stochastic process model of airport pavement performance degradation was estab- 
lished. The transition probability of the semi-Markov process model was estimated by using the Weibull 
function parameter model of survival analysis. In the process of building the model ,the dual effects of flight 
traffic and pavement surface thickness were considered ,and the measured pavement performance data set 
of local civil airports was adopted to analyze the influence of these two factors in Henan Province from 
2007 to 2017. The performance degradation process of runways of the civil airport in Henan province was 
predicted by using the transition probability of the semi-Markov process model. The transition probability 
was updated by adopting the Bayesian statistical analysis method based on Markov Chain Monte Carlo 
(MCMC) simulation technology with Metropolis-Hasting ( MH) sampling method ,using a priori informa- 
tion obtained from the irregular detection data of the airport runway. Finally ,the updated model is applied 
to predict the performance of the local civil airport runway in Henan Province. The results show that the 
Bayesian analysis method based on MCMC can fully integrate the regular and irregular detection data to 
determine the prior information of parameters of the dynamic semi-Markov process model,and can effec- 


=stively update the performance prediction model of the airport runway and improve the prediction accuracy 


CNof the model. 
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CO 机场 道 面 性 能 状况 一 直 都 是 机 场 设施 运 维 部 门 
最 奖 注 的 问题 ,通过 道 面 性 能 预测 ,可 以 对 养护 维修 
决策 提供 指导 ,有 助 于 将 道 面 性 能 维持 在 适当 的 水 
琵 芒 期 以 来 道 面 性 能 预测 都 是 机 场 道 面 领域 研究 
的 重点 。 

.二 自 2002 年 以 来 ,同济 大 学 结合 上 海 机 场 道 面 调 
查 评 价 与 管理 的 实践 ,在 借鉴 国外 机 场 道 面 管理 系 
统 的 基础 上 , 开发 建立 了 上 海 机 场 道 面 管理 系统 
(SEAPMS) 1 。 目 前 ,国内 的 大 型 民航 枢纽 机 场 已 
经 俩 用 机 场 道 面 管理 系统 对 机 场 士 木 工程 设施 进行 
检测 ,评价 以 及 运营 维护 1。 总 体 来 说 ,国内 道 面 管 
理 系统 在 国内 机 场 跑道 养护 、 维 修 和 管理 中 的 应 用 
越 来 越 广 泛 。 通 过 近 20 年 的 机 场 道 面 管理 系统 的 
应 用 ,中 国民 航 管理 部 门 也 积累 了 较 多 的 数据 资 
料 55 ,包括 每 5 年 对 机 场 道 面 进行 定期 检测 的 数 
据 日常 运 维 数据 以 及 道 面 修缮 和 扩建 数据 等 。 目 
前 虽然 有 研究 人 员 :“1 采 用 这 些 数据 进行 分 析 和 统 
计 , 并 对 机 场 道 面 管理 提供 了 一 定 的 技术 指导 ,但 是 
相 比 于 国内 的 路 面 管理 系统 的 数据 库 来 说 ,机 场 道 
面 管理 系统 数据 库 的 利用 率 仍然 处 于 较 低 水 平 。 
国外 应 用 路 面 管理 系统 较 早 ,建立 了 适用 各 自 
管理 系统 的 模型 ,包括 线性 外 推 模型 “S” 型 曲线 模 
型 .约束 多 项 式 模 型 等。 而 应 用 于 机 场 道 面 使 用 
性 能 预测 中 最 具 代 表 性 的 经 验 模型 为 PAVER 道 面 


管理 系统 模型 ,Shahin 等 ”在 道 面 状况 指数 (pave- 
ment condition index,PCI) 基础 上 建立 和 发 展 了 基于 
“ 道 面 族 "方法 的 回归 模型 ,该 模型 以 PCI 为 预测 指 
标 ,以 多 项 式 作为 函数 形式 ,通过 最 小 二 乘法 进行 参 
数 估 计 。 

相 比 经 验 模型 而 言 ,力学 模型 ”具有 明确 的 理 
论 基础 ,为 了 兼顾 力学 和 经 验方 法 的 优点 ,近年 来 ， 
越 来 越 多 的 学 者 倾向 于 力学 -经 验 模型 。1995 年 
FAA 提出 以 累积 疲劳 损伤 原理 、 弹 性 层 状 体系 理论 
和 通行 -覆盖 率 为 基础 ,针对 新 一 代 大 型 飞机 的 道 面 
设计 方法 ,并 以 该 设计 方法 的 关 过 程 作为 道 面 结构 
剩余 寿命 的 预 估 模 型 ,该 模型 是 机 场 道 面 力学 - 
经 验 模 型 的 典型 代表 。 同 济 大 学 借鉴 设计 过 程 逆转 
的 思想 ,在 我 国 机 场 道 面 设计 规范 的 基础 上 ,建立 了 
我 国 机 场 刚性 道 面 结构 剩余 寿命 的 预 估 模 型 ,并 
在 国内 多 个 民用 机 场 的 结构 评价 中 应 用 和 验证 。 

道 面 性 状 的 预测 有 很 大 的 不 确定 性 ,确定 型 模 
型 无 法 解决 此 问题 ,所 以 目前 不 少 铺面 管理 系统 应 
用 概率 方法 开发 了 道 面 性 状 预测 模型 。 概 率 型 预测 
模型 包括 残存 曲线 概率 模型 、 贝 叶 斯 概率 模 
型 半 马 尔 可 夫 概 率 模型 ”和 马尔 可 夫 概 率 模 
型 ,其 中 应 用 最 广泛 、 最 完善 的 是 马尔 可 夫 模 型 。 

由 于 机 场 道 面 实施 监测 的 成 本 较 高 ,目前 机 场 
道 面 系统 数据 库 中 检测 数据 占 比 较 大 。 这 些 检 测 数 
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据 中 隐藏 着 大 量 信 息 , 因 此 有 必要 对 这 些 信息 进行 
数据 挖掘 ,并 融合 新 检测 数据 对 原 有 数据 进行 分 析 ， 
获取 有 用 信息 ,以 指导 道 面 运 维和 管理 。 因 此 ,本 研 
究 以 国内 机 场 道 面 性 能 评 佑 方法 为 基础 ,从 某 民 航 
机 场 定 期 检测 数据 中 提取 样本 ,建立 机 场 道 面 性 能 
预测 模型 ,为 机 场 跑 道 养 护 决策 提供 技术 支持 ,同时 
也 为 机 场 管理 部 门 合理 地 分 配 运 维 资金 提供 参考 。 


1 机 场 跑 道 性 能 退化 模型 的 建立 


1.1 半 马 尔 可 夫 过 程 模型 
在 机 场 道 面 性 能 退化 的 过 程 中 , 如 果 机 场 道 面 
区 域 的 未 来 状况 依赖 于 道 面 区 域 的 当前 状况 , 则 该 
随机 过 程 满足 马尔 可 夫 性 。 由 于 道 面 性 能 衰减 过 程 
中 状态 和 时 间 都 是 连续 变量 ,因此 ,在 应 用 马尔 可 夫 
预测 模型 时 常 做 一 些 规定 ”。 机 场 道 面 退化 过 程 
其 闫 态 持续 时 间 具 有 明显 的 半 马 尔 可 夫 性 "。 因 
由 森 研 究 利 用 半 马 尔 可 夫 模 型 分 析 机 场 道 面 退 化 
的 风机 过 程 。 
CD 道 面 性 能 退化 的 半 马 尔 可 夫 预 测 模型 建立 步 又 
如 肌 。 
CNI ) 选择 道 面 状况 指数 PCI(0 ~ 100 ) 作为 预测 指 
梳 S 刘 分 指标 状态 空间 。 将 道 面 PCI 离散 为 由 优 、 
良 -中 次 , 差 5 个 状态 组 成 的 状态 空间 。 设 任意 时 
刻 匡 道 面 性 能 状况 处 于 各 种 状态 的 概率 分 布 为 
XK ,概率 转移 的 观测 周期 以 年 为 单位 。 过 程 的 初 
始 醚 刻 记 为 1, 则 道 面 性 能 的 初始 分 布 向 量 为 
Xo) = (Qo ,ao ，… ,ao5 ) ， > ao =1。6 


2) 建立 半 马 尔 可 夫 转 移 概率 矩阵 。 半 马尔 可 夫 
链 主 要 通过 转移 概率 矩阵 预测 道 面 性 能 随时 间 的 训 
变 过 程 。 假 设 在 日 常 维修 养护 条 件 下 , 道 面 使 用 性 
能 由 低 水 平 状态 向 高 水 平 状态 转移 的 情况 不 会 发 
生 。 即 当 ;i >7 时 六 =0; 则 转移 矩阵 可 表示 为 
[pn Pi P13 Pa pis | 


0 po Pp» Pa Ps 
P=| 0 0 ps Pa Pas (1) 
0 0 0 py pas 
10 0 0 0 ps 
3) 在 实际 预测 过 程 中 ,考虑 到 半 马 尔 可 夫 链 的 
无 后 效 性 , 由 Chapman-Kolmogorov 等 式 可 得 到 由 
时 刻 状态 到 任意 上 时 刻 状态 的 转移 矩阵 己 ( 切 如 式 
(2) 所 示 。 根 据 初始 状态 的 概率 向 量 导 (0) , 则 半 马 


广 hinsyiv 全 人 作 革 日生 | 
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尔 可 夫 预 测 模型 可 写 为 X(t) = 四 (0)PCi) 。 


[Pu(t) polt) pa pa pis(t)| 

0 p(t) pa(t) pult) pzs(t) 

P(i)=| 0 0 pa(t) palt) pas(t) 
0 0 0 palt) pas(i) 

| 0 0 0 0 pss(t) | 


(2) 

1.2 转移 概率 的 计算 方法 

对 道 面 性 能 退化 过 程 的 参数 模型 进行 估计 ,可 
以 计算 半 马 尔 可 夫 过 程 的 转移 概率 。 一 般 情况 下 ， 
离散 状态 退化 模型 的 特征 是 : 某 一 状态 的 转移 概率 
在 一 个 检查 周期 A(1 年 或 2 年 ) 内 为 保持 不 变 。 在 
时 刻 t 道 面 性 能 状态 等 级 变量 Y=1 的 条 件 下 ,下 一 
时 刻 #+A4 观测 到 道 面 性 能 状态 等 级 变量 Y=0 的 
概率 是 条 件 概率 ,用 R(i,A) 表 示 , 可 以 由 式 (3) 计 
算 , 式 (3) 中 Probp(，) 表 示 概 率 ,了 表示 时 间 变 量 。 

R(t1,A) =Prob(1<T<t+AIT>t) 


_ Prob(1 <T <i+A) 
Prob(T >1) 03) 
R(1,A) 可 以 用 生存 分 析 "” 分 布 函 数 (1) = 
P(Ts<) 或 生存 函数 S(t) = PC(7 > 发 计算 。 


依据 生存 分 析 理 论 可 以 得 到 


SG) =1-P(D) =1- [f(a (4) 


式 中 f(t) 为 随机 变量 1 的 概率 密度 函数 。 

可 以 使 用 各 种 不 同 的 概率 分 布 来 描述 失效 概 
率 , 威 布尔 分 布 孔 数 比 较 适合 描述 土木 工程 结构 的 
失效 概率 和 。 因 此 ,本 研究 利用 威 布尔 分 布 概率 密 
度 函 数 计算 生存 函数 5S(1)。 威 布尔 分 布 概率 密度 
函数 可 以 定义 为 

fliilp,A) =Ap(At)’ exp( ~- (At)’) (5) 

式 中 pW 为 威 布尔 分 布 函 数 的 参数 。 

t 时 刻 的 瞬时 失效 概率 可 以 用 危险 率 函 数 入 (1) 
表示 , 即 


A(i) =Ap (AL)”! =p 和 PiP-1 (6) 


机 场 跑道 的 性 能 退化 是 多 种 因素 (比如 设计 特 
性 、 飞 机 交通 荷载 .使 用 环境 、 使 用 年 限 和 维修 历史 
等 ) 综 合作 用 的 结果 ,因此 可 以 用 一 组 影响 因素 (或 
解释 变量 ) 对 跑道 的 性 能 退化 过 程 进行 解释 。 通 过 
对 式 (6) 引 入 解释 变量 ,可 以 将 其 转化 为 包含 待 佑 
参数 的 形式 ,为 了 保证 危险 率 函数 的 非 负 性 ,采用 指 


数 函 数 引 入 解释 变量 , 即 
A=e Bs (了 7) 
式 中 :于 为 外 生 解 释 变 量 列 向 量 (其 中 常数 项 对 应 分 
量 值 1) ;有 为 待 估 参 数 行 向 量 。 
将 式 (7) 代 入 式 (6) ,并 进行 积分 可 以 得 


A(1) = (AL)? (8) 
进而 可 以 得 到 生存 函数 S(t) 的 表达 式 , 即 
S(ti) =exp ( - (Ai)’) (9) 


根据 式 (3) 和 式 (9) 可 以 得 到 在 任何 时 刻 1, 经 
过 一 段 时 间 A 从 状态 i 转移 到 低 一 级 状态 〖 的 概率 
pa(i,A), 即 


P (LA) =R(,A) _Prob(t<T <t+A) 


Prob(T >1) 
_Fi(t+A) -F(t) _Fi(t+A) -F(t) 
1-F(t) Si(1) 


CI (10) 
式 趾 :F(1) 为 在 状态 i 的 持续 时 间 变量 1 的 累积 分 
布 亢 数 ;Si;(1) 为 状态 i 持续 时 间 的 生存 函数 。 

一 对 于 半 马 尔 可 夫 过 程 中 的 转移 概率 ,本 研究 基 
于 文献 [20] 提出 的 方法 来 分 析 多 状态 情况 下 ( 优 、 
让 Gh 次 、 差 ) 机 场 道 面 性 能 退化 过 程 的 转移 概率 。 
CN 考虑 到 机 场 道 面 性 能 退化 过 程 为 半 马 尔 可 夫 随 
机 时 程 ,在 +=0 时 刻 之 前 , 道 面 在 状态 i 持续 的 时 间 
用 心 表示 ,其 下 一 个 可 能 转移 的 目标 状态 为 状态 让 
准 寻 系统 历经 中 间 状 态 k。 其 中 时 间 5 为 道 面 在 状 
态 c6 条 件 下 , 道 面 的 使 用 年 数 (用 Age 表示 ) 。 
.和 对 于 半 马 尔 可 夫 链 , 其 转移 概率 方程 为 

SD) i=] 

py 六 (CORP —%) ,i zj 


式 中 :为 状态 iV 之 间 的 中 间 状 态 ;fi(%) 为 时 间 
开始 ,x 时 刻 发 生 第 一 次 转移 时 状态 i 的 持续 时 间 概 
率 密度 函数 ;Py(1-*) 是 在 剩余 时 间 (! -*) 内 从 状 
态 上 转移 到 状态 j 的 概率 ,可 利用 式 (10) 计 算 。 
在 系统 进入 状态 的 条 件 下 ,在 下 一 个 时 间 间 
隔 A=t-x 内 ,状态 上 将 要 转移 为 状态 j 的 概率 可 以 
用 状态 j 转移 发 生 之 前 的 状态 的 持续 时 间 黑 积分 
布 函数 F(t) 计算。 将 1=0,A =1-« 代 入 式 (11)， 
可 得 到 
| Fy(t ~ «) ~ Fy(0) 
.I 
j=k+l (12) 
1.3 ”概率 型 退化 过 程 的 更 新 方法 


© P00) - (11) 
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为 了 解决 传统 回归 模型 更 新 方法 的 缺点 , 目前 
发 展 了 一 种 更 新 模型 参数 的 贝 叶 斯 ( Bayesian ) 统计 
方法 。 近 年 来 ,许多 研究 者 应 用 贝 叶 斯 方法 对 基础 
设施 性 能 退化 模型 的 参数 进行 了 更 新 ” 。 

因为 机 场 道 面 性 能 退化 过 程 的 危险 率 随 使 用 时 
间 的 增长 逐渐 增加 ,所 以 机 场 跑道 的 性 能 退化 可 以 
采用 式 (7) 表 示 的 危险 率 函 数 和 (i) ,并 采用 生存 分 
析 的 参数 模型 方法 求 出 危险 率 函数 和 (i) 的 参数 方 
程 ,利用 1.2 节 所 述 的 方法 求 出 任意 时 刻 1 和 任意 
周期 A( 本 研究 采用 机 场 道 面 定 期 检测 周期 5 年 ) 的 
状态 转移 概率 数值 ,然后 根据 式 (2) 预测 道 面 使 用 
性 能 任意 时 刻 各 状态 的 分 布 概 率 。 

生存 分 析 利 用 机 场 道 面 定期 检测 数据 样本 组 成 
的 历史 数据 集 进行 分 析 , 得 出 各 状态 的 转移 概率 ,对 
此 模型 进行 分 析 评 价 。 

然后 ,将 新 收集 的 不 定期 检测 数据 添加 到 历史 
数据 中 ,并 使 用 MCMC'” 对 生存 分 析 参 数 模型 进行 
仿真 并 估计 结果 。 当 模型 结果 更 新 后 ,采用 贝 叶 斯 
生存 分 析 参 数 模型 得 到 的 结果 预测 机 场 道 面 的 性 能 
退化 。 最 后 ,将 生存 分 析 参 数 模 型 的 结果 与 基于 
MCMC 仿真 的 后 验 估计 结果 进行 比较 ,从 而 评价 由 
叶 斯 更 新 方法 的 效果 。 


2 机 场 道 面 数 据 的 准备 


机 场 道 面 历史 数据 是 建立 机 场 道 面 性 能 退化 预 
估 模 型 的 基础 ,近年 来 ,我 国 已 积累 了 近 60 多 个 民 
用 机 场 的 道 面 PCI 测试 数据 ,这 些 数据 一 方面 能 够 
代表 我 国民 用 机 场 道 面 的 性 能 状况 , 男 一 方面 也 能 
代表 道 面 性 能 现场 测试 数据 的 典型 特征 。 本 研究 以 
两 组 不 同 特点 的 机 场 道 面 测试 历史 数据 为 数据 源 ， 
数据 集 1 是 我 国民 用 机 场 道 面 PCI 现场 测试 历史 数 
据 集 , 数 据 集 2 是 某 民 航 机 场 跑道 2007 年 ~2017 年 
现场 测试 数据 集 。 其 中 数据 集 1 的 特点 和 详细 信息 
参见 文献 [22]。 这 里 只 对 数据 集 2 进行 简要 介绍 。 

数据 集 2 是 某 民 航 机 场 2007 年 .2013 年 和 
2017 年 南 飞行 区 跑道 滑行 道道 面 的 检测 时 间 序 列 
数据 。 该 机 场 于 1997 年 建成 通航 ,至 今 其 南 飞 行 区 
道 面 已 使 用 23 年 ,已 接近 设计 使 用 寿命 。 为 了 全 面 
掌握 南 飞 行 区 道 面 综 合 情 况 ， 机 场 管 理 部 门 于 
2007 年 、2013 年 和 2017 年 分 别 对 机 场 南 飞行 区 跑 
道 、 滑 行道 和 联络 道道 面 进行 了 全 面 检测 ， 并 对 现 
场 检测 数据 进行 了 综合 分 析 ， 形 成 了 具有 较 高 价值 
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的 数据 库 。 
表 1 机 场 道 面 历史 数据 内 容 简 表 


Tab. 1 Summary table of airport pavement history data set content 


ee 

来 源 持续 时 间 /年 数 /年 交通 量 em 
1 1 10 1 8 35 934 34 
和 2 1 10 0 8 35 934 32 
3 1 12 0 12 5 132 2 
4 1 12 1 10 71 236 38 
20 2 12 0 12 74 325 38 
21 2 10 1 8 74 325 38 
33 2 14 1 12 74 325 34 
34 2 12 0 12 74 325 34 


三 本 研究 采用 《民用 机 场 道 面 评价 技术 管理 规 


薄 兴 MHZT5024 - 2019 ) 的 方法 以 道 面 评价 区 域 的 
“ 道 面 状况 指数 PCI 平均 值 "来 评定 道 面 表 观 状况 等 
级 .具体 等 级 划分 :PCI > 85 为 优 ;70 <PCI < 85 为 
旗 335 <PCI <70 为 中 ;40<PCI < 55 为 次 ;PCI <40 
洲 鸡 .数据 中 , 道 面 个 体 样本 所 处 等 级 为 优 . 良 、 中 、 
次 而 种 状态 (其 差 状 态 等 级 为 退化 过 程 的 结束 状 
态 六 因此 本 研究 主要 对 这 4 种 状态 进行 分 析 , 估 计 
持 肥 时 间 模型 并 计算 其 转移 概率 。 

“在 生存 分 析 中 , 感 兴趣 的 因 变量 是 道 面 在 给 定 
条 传 下 维持 某 种 等 级 状态 的 时 间 。 这 个 变量 定义 为 
状态 时 间 (time-in-state,TIS) 。 对 于 TIS 观测 值 为 右 
删 路 ”的 道 面 ,由 于 其 进入 当前 状态 等 级 的 时 间 未 
知 5 只 知道 发 生 在 检测 开始 时 间 之 前 的 某 个 时 间 。 
六 方便 分 析 , 对 于 右 删 失 的 道 面 个 体 的 使 用 年 数 
采用 检测 时 刻 的 年 数 近似 作为 进入 当前 状态 等 级 的 
使 用 年 数 。 按 照 上 述 方法 ,最 后 得 到 107 组 机 场 首 
面 的 生存 分 析 数据 ,如 表 1 所 示 ( 只 显示 部 分 样本 ) 。 
其 中 第 1 ~ 16 组 样本 为 数据 集 1 ,涵盖 的 机 场 刚性 道 
面 厚度 范围 为 18 cm ~ 40 em, 基本 黎 盖 了 我 国 不 同 
等 级 民用 机 场 道 面 结构 设计 厚度 ,飞行 交通 量 范围 
为 100 ~71 236 ,使 用 年 限 范围 为 15 年 ~23 年 。 第 
17 ~ 107 组 样本 为 数据 集 2 ,数据 集 2 是 某 民 航 机 场 
2007 年 .2013 年 和 2017 年 南 飞行 区 跑道 滑行 道道 


| 
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次 道 面 大 修 或 者 建成 投入 使 用 时 为 起 点 直至 退化 到 
当前 状态 的 时 间 。 飞 行 交 通 量 为 当前 状态 转移 时 机 
场 运行 的 主要 机 型 的 累积 换算 架次 。 厚 度 为 道 面 水 
泥 混 凝 土 面 层 厚度 ,单位 为 cm。 


3 结果 分 析 


3.1 半 马 尔 可 夫 过 程 模型 结果 

根据 文献 [23] 以 及 前 文 的 分 析 , 具体 的 影响 因 
素 考 虑 飞行 交通 量 \ 道 面 厚度 和 使 用 年 数 三 者 ,其 中 
飞行 交通 量 采 用 当前 状态 转移 时 机 场 运 行 的 主要 机 
型 的 累积 换算 架次 的 对 数值 。 危 险 率 函数 表达 为 

a i 
nA =B +B,lgN +BsH +BaAge 
式 中 :B, 为 影响 因素 的 待 估计 系数 ;Age 为 使 用 年 数 ， 
是 从 最 近 一 次 道 面 大 修 或 者 建成 投入 使 用 时 为 起 点 
直至 退化 到 当前 状态 的 时 间 ;lgN 为 飞行 交通 量 的 对 
数值 ;五 为 机 场 道 面 面 层 混凝土 厚度 ,单位 为 cm。 

4 种 状态 (优良 .中 次 ) 等 级 的 持续 时 间 函 数 
参数 佑 计 的 结果 如 表 2 所 示 。 表 2 中 的 系数 估算 值 
均 在 显著 性 水 平 0.05 上 具有 统计 学 意义 ;“* * ” 
表示 系数 估算 值 没 有 显著 的 统计 学 意义 ; 插 号 内 的 
数值 为 显著 性 水 平 0. 1 的 边际 显著 卡 方 统计 值 ;p 
为 威 布尔 分 布 函 数 的 形状 参数 。 

从 表 2 参数 估计 可 以 得 出 ,除了 状态 等 级 2 
( 良 ) 的 使 用 年 数 Age 具有 显著 统计 意义 外 ,其 他 状 
态 等 级 的 使 用 年 数 Age 均 没 有 显著 的 统计 意义 ,其 
他 所 有 影响 因素 的 系数 估计 值 的 显著 性 水 平均 小 于 
0.1 ,这 说 明 飞 行 交 通 量 和 道 面 厚度 对 于 模型 具有 显 
著作 用 。 状 态 等 级 1( 优 ) 和 状态 等 级 3( 中 ) 的 飞行 
交通 量 lgN 的 系数 估计 值 为 负 值 , 即 飞 行 交 通 量 和 
危险 率 为 负 相 关 , 说 明 随 着 飞行 交通 量 的 增加 ,人 危险 
率 逐 渐 增 加 。 换 言 之 , 当 飞 机 的 起 降 架 次 增加 时 , 飞 
机 和 荷载 作用 在 机 场 道 面 的 频率 增加 ,在 飞机 人 荷载 的 
重复 作用 下 ,机 场 跑道 结构 的 应 力 会 出 现 累 积 , 当 达 
到 一 定数 值 后 ,会 引起 机 场 道 面 板 在 应 力 薄弱 部 位 
产生 裂缝 ,进一步 引起 PCI 数值 的 减 小 , 则 机 场 道 面 
维持 当前 状态 等 级 的 概率 减 小 ,从 而 机 场 道 面 维持 


面 的 检测 数据 ,刚性 道 面 厚度 为 34 cm( 跑道 中 部 ) 、 
38 cm (跑道 两 端 和 滑行 道 ) ,飞行 交通 量 为 74 325， 
使 用 年 限 为 19 年 。 

表 1 中 删 失 与 耕 ( censor) ,其 中 “1" 表示 右 册 
失 ,“0" 表 示 状 态 转移 发 生 。 使 用 年 数 是 从 最 近 一 


当前 状态 等 级 的 持续 时 间 TIS 的 危险 增加 , 其 转移 
到 下 一 状态 等 级 的 概率 增 大 。 

男 外 ,状态 等 级 优 和 中 的 道 面 厚度 五 的 系数 估 
计 值 为 正 值 , 即 道 面 厚度 和 和 危险 率 为 正 相 关 , 说 明 随 
着 道 面 厚度 的 增 大 ,人 危险 率 逐 渐 减 小 。 换 言 之 , 随 着 


机 场 道 面 厚度 的 增 大 ,机 场 跑道 结构 的 刚度 增 大 ,其 
结构 性 承载 力 增强 ,能 够 承担 的 飞机 荷载 作用 等 级 
以 及 作用 频次 也 相应 增加 ,机场 道 面 维持 当前 状态 
等 级 的 概率 增 大 ,从 而 维持 当前 状态 等 级 的 持续 时 
间 TIS 的 危险 越 小 ,其 转移 到 下 一 状态 等 级 的 概率 
也 越 小 。 

对 于 状态 等 级 良 的 系数 估计 值 ,除了 截 距 系数 
B1 和 道 面 厚度 及 的 斜率 系数 B; 显 著 性 水 平 为 边际 显 
著 ( 显 著 水 平 0.10) 外 ,其 他 影响 因素 的 系数 估计 值 
的 显著 性 水 平均 小 于 0.05, 这 说 明 模 型 考虑 的 变量 
具有 显著 的 作用 。 由 于 飞行 交通 量 lgN 的 估计 系数 
为 正 值 ,与 状态 1 的 估计 值 ( 负 值 ) 正 负 性 相反 ,这 并 
不 能 说 明 对 于 机 场 道 面 状态 等 级 2 的 性 能 退化 过 程 
和 状态 等 级 1 的 退化 过 程 的 影响 因素 不 同 ,而 可 能 
是 由 于 其 他 解释 变量 ( 比如 使 用 年 数 ) 的 内 在 影响 
作用 掩盖 了 飞行 交通 量 leN 的 作用 。 
< 二 此 外 ,对 于 状态 等 级 次 的 系数 估计 值 ,飞行 交通 
量 i8N 的 系数 B, 为 正 值 , 道 面 厚度 及 的 系数 Bs 为 负 
值 ,这 和 状态 等 级 1 或 3 的 相应 系数 正 负 正好 相反 。 
这 些 现象 并 不 能 说 明 飞 行 交通 量 和 道 面 厚度 对 于 状 
巷 尝 级 4( 次 等 ) 持 续 时 间 的 影响 和 其 他 的 状态 等 级 
是 相反 作用 。 分 析 本 研究 所 采用 的 数据 集 可 以 发 
更; 机 场 道 面 在 退化 状态 等 级 为 次 等 时 ,其 PCI 数值 
小 各 55 ,这 说 明 机 场 道 面 的 功能 性 表 观 状况 已 经 影 
响 到 机 场 道 面 的 内 在 结构 物理 性 能 (应 力 、 应 变 
等 呵 此 时 机 场 道 面 性 能 退化 过 程 不 只 与 飞行 交通 
量 S 道 面 厚 度 有 关 , 还 与 其 他 影响 因素 ( 比如 环境 温 
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度 湿度 . 道 面 材料 .跑道 结构 组 成 等 ) 有 关系 。 至 于 
这 些 影 响 因 素 的 作用 规律 如 何 , 是 下 一 步 需要 研究 
的 内 容 。 

道 面 退化 过 程 的 初始 时 刻 记 为 为 , 道 面 性 能 的 
初始 分 布 向 量 为 Cc" (10) =〈(aoy，ao，…，ao5 ) ， 
> ,a0 =1 ;其 中 au 是 机 场 道 面 在 加 时 刻 处 于 状态 
i 的 概率 ,i =1,2,…,5。 由 Chapman-Kolmogorov 等 
式 可 得 到 任意 上 时刻 状 态 的 转移 矩阵 书 等 于 初始 分 
布 向 量 C" (io) 与 预测 转移 概率 和 矩阵 P;(1) 的 乘积 即 
P=C'(t) =C (to)P;(it)。 根据 《民用 机 场 道 面 评 
价 技术 管理 规范 》 的 状态 评价 等 级 PCI 数值 向 量 
Po = (100 ~85,85 ~70,70 ~55,55 ~40,40 ~0) 可 
以 得 出 1 时 间 后 机 场 道 面 的 PCI 预测 值 为 Po(t) = 
C(t) xPuo 

取 初 始 状态 向 量 C(t。) = (10 000) ,分 别 以 数 
据 集中 第 8 组 数据 和 第 22 组 数据 作为 分 析 案 例 , 应 
用 半 马 尔 可 夫 过 程 模 型 得 到 各 年 的 状态 转移 概率 算 
阵 ,计算 得 到 该 组 道 面 的 使 用 性 能 概率 预测 结果 。 
其 中 第 8 组 样本 的 转移 概率 和 矩阵 图 见 图 1 和 图 2 所 
示 ( 其 中 “P;” 表 示 第 8 组 样本 的 模型 状态 i 转移 到 
状态 j 的 概率 )。 两 组 道 面 的 使 用 性 能 预测 结果 如 
图 3 所 示 。 图 3 中 , 蓝 色 实 线 代表 机 场 道 面 PCI 值 
状态 区 间 中 值 预 测 曲 线 , 上 下 虚线 分 别 为 PCI 上 限 
值 和 下 限 值 预测 曲线 , 带 三 角 标 志 的 点 划 线 为 机 场 
道 面 PCI 实测 值 ”DatalS08 "表示 数据 集 1 中 的 第 8 
组 样本 ,其 他 编号 类 似 。 


表 2 状态 持续 时 间 威 布尔 函数 参数 估计 值 


Tab.2 Parameter estimate of Weibull function of state duration 


优 (1) 良 (2) 中 (3) 次 (4) 
系数 EE- 一 下 一 » 一 Wi 到 .一 > 
估算 值 标准 误差 估算 值 标准 误差 估算 值 标准 误差 估算 值 标准 误差 
Bi 2.93 3.25x10-1 2.96x10-1(2.65) 1.82x10-! 1.03 1.91 x10-1 3.26 2.66 x10-! 
Bb, -1.83x10-1 8.28 x107? 2.38 x10-1 5.16xl10-” -1.75x10-! 5.90x10-* 3.90x10-! 9.97x10™ 
Bs 2.21xl0-” 1.08x10-? 1.63x10-2(3.04) 9.30x10-3 9.59x10-2 9.10x10-3 -6.96x10”? 1.25 x10-? 
Ba 六 六 六 六 2.96 x10-? 7.90 x10-3 六 六 六 六 六 六 六 六 
1/p 1.18 x10-! 1.86x10-? 1.15 x10-1 1.18x10-” 1.23x10-! 1.32x10-* 2.19x107! 2.24 x107™ 
p 8.50 1.34 8.72 8.94 x10-! 8.15 8.79 x107! 4.57 4.69 x10-1 


由 图 3 可 知 ,实测 道 面 PCI 数据 基本 均 分 布 在 
模型 预测 曲线 的 预测 区 间 之 内 ,这 说 明 半 马尔 可 夫 
过 程 模型 对 于 两 组 (第 8 组 和 第 22 组 ) 道 面 PCI 数 
据 均 能 较 好 的 拟 合 。 另 外 从 图 3(a) 和 图 3(p) 的 比 
较 可 以 得 出 : 半 马 尔 可 夫 模 型 对 于 机 场 道 面 性 能 退 


化 过 程 的 预测 曲线 ,在 退化 过 程 前 期 ,预测 区 间 的 上 
限 接近 实测 值 ,而 在 退化 过 程 中 期 和 后 期 ,预测 区 间 
的 中 值 接近 实测 值 ; 即 半 马 尔 可 夫 过 程 模型 对 于 机 场 
道 面 退化 过 程 的 状态 等 级 3 ~5( 中 ~ 差 ) 的 持续 时 间 
预测 效果 要 好 于 状态 等 级 1( 优 ) 或 2( 良 )。 主 要 原因 
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可 能 和 道 面 性 能 退化 的 力学 特点 有 关 , 在 道 面 使 用 的 
前 期 ,由 于 使 用 时 间 较 短 ,飞机 荷载 作用 次 数 较 少 , 道 
面 的 结构 性 承载 能 力 基本 维持 在 设计 水 平 ,保持 在 优 
良 等 级 的 概率 比较 大 ,从 而 其 PCI 数值 基本 维持 在 高 
位 水 平 。 但 是 在 道 面 使 用 的 中 后 期 ,由 于 使 用 时 间 的 
增加 ,飞机 荷载 作用 次 数 的 累积 , 道 面 的 结构 性 承载 
能 力 逐 渐 被 消耗 , 道 面 使 用 状态 等 级 转移 的 概率 增 
大 ,从 而 其 PCI 数值 更 接近 于 中 低位 水 平 。 

从 状态 等 级 持续 时 间 上 分 析 , 半 马尔 可 夫 过 程 
模型 对 于 机 场 道 面 PCI 退化 到 50( 第 8 组 道 面 ) 或 
退化 到 70( 第 22 组 道 面 ) 的 持续 时 间 的 预测 值 均 较 
实际 值 大 ( 预测 曲线 均 在 实测 值 曲线 的 右 侧 ) ,这 说 
明 半 马 尔 可 夫 过 程 模型 对 于 机 场 道 面 退化 过 程 的 巴 
测 相对 保守 ,这 和 半 马 尔 可 夫 过 程 模型 考虑 了 飞行 
交通 量 和 道 面 厚度 的 影响 因素 、 从 而 使 预测 数据 能 
较 妈 地 反映 实际 的 道 面 性 能 退化 过 程 有 一 定 关系 。 
3 史 - 机 场 道 面 性 能 退化 模型 的 更 新 结果 
写 如 果 马 尔 可 夫 链 具有 遍历 性 且 数 据 链 足够 长 ， 
则 得 估 参 数 后 验 分 布 将 会 接近 于 其 真实 分 布 。 一 般 
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情况 下 , 待 估 参 数 的 先 验 分 布 可 以 任意 选择 。 然 而 ， 
为 了 获得 贝 叶 斯 分 析 的 良好 收敛 性 能 , 则 要 求 先 验 
分 布 的 选择 形式 应 接近 于 参数 的 后 验 分 布 。 
此 ,本 节 采 用 的 MCMC 的 贝 叶 斯 分 析 模 型 ,根据 3.1 
节 中 生存 分 析 参 数 模 型 的 估计 结果 表 2 设置 先 验 分 
布 。 此 外 ,本 研究 假设 等 估 的 先 验 分 布 是 单 变量 正 

为 了 验证 本 研究 所 提出 的 MCMC 贝 叶 斯 模型 
的 精度 和 有 效 性 ,将 新 收集 的 91 组 机 场 道 面 性 能 类 
据 样 本 加 入 到 历史 数据 中 。 根 据 表 2 的 参数 估计 结 
果 , 确 定式 (13) 待 佑 参数 B; 的 先 验 分 布 为 正 态 概率 
分 布 ,形状 参数 p 的 先 验 分 布 为 伽 马 分 布 。 

由 于 机 场 道 面 性 能 退化 过 程 是 五 状态 空间 的 半 
马尔 可 夫 链 ,所 以 采用 贝 叶 斯 理论 更 新 模型 针对 5 
个 状态 ( 优 、 良 .中 次 、 差 ) 分 别 分 析 。 根 据 表 2 中 数 
据 和 新 增加 的 机 场 道 面 性 能 检测 数据 ,确定 各 个 参 
数 的 先 验 分 布 ,如 表 3 所 示 。 表 中 N(3 ,0.3 ) 表 示 均 
值 为 3、 标 准 差 为 0.3 的 正 态 分 布 ;G(0. 02,0. 01) 
表示 尺度 参数 0. 02 ,形状 参数 0.01 的 伽 玛 分 布 。 


表 3 MCMC 方法 贝 叶 斯 更 新 后 的 模型 参数 值 
Tab.3 Model parameter values updated by using Bayesian approach with MCMC method 


CY 马尔 可 大 95% 置信 区 间 
系数 先 验 分 布 估计 值 标准 差 
CN 标准 误差 下 限 上 限 
状态 1( 优 ) 
GB N(3,0.3) 2.96 0.19 0.005 40 2.62 3.37 
> N( -0.2,0.08 ) 一 0.17 0.048 0.00191 一 0.27 一 0.07 
FF N(0.01,0.01) 0.017 0.006 0.000 22 0.004 0.029 
人 C(8,0.8) 8.64 1.18 0.0391 6.48 11.16 
状态 2( 良 ) 
Bi N(0.3,0.2) 0.32 0. 133 0.004 2 0.068 0.59 
Bb N(0.2,0.05) 0.22 0.035 0.0015 0.15 0.29 
Bs N(0.02,0.01) 0.02 0.0059 0.000 28 0.009 0.031 
Ba N(0.02,0.01) 0.027 0.0058 0.000 18 0.016 0.038 
p G(8,0.8) 8.59 0.85 0.028 6.91 10.22 
状态 3( 中 ) 
Bi N(0.5,0.1) 0.62 0.078 0.0024 0.47 0.77 
Bs N( -0.2,0.06) 一 0.17 0.038 0.0016 一 0.24 一 0.096 
Bs N(0.1,0.001) 0.099 0.0055 0.000 23 0.088 0.11 
p G(8,2) 6.76 0.67 0.02 5.43 8.09 
状态 4( 次 ) 
Bi N(3,0.3) 3.36 0.185 0.0058 3.00 3.73 
Bs N(0.4,0.08) 0.45 0.047 0.0016 0.36 0.54 
Bs N( -0.1,0.01) 一 0.083 0.0069 0.000 22 一 0.097 一 0. 069 
p G(4,1) 4.28 0.43 0.0145 3.48 5.17 


根据 上 述 待 估 参 数 的 先 验 分 布 ,采用 MCMC 的 
贝 叶 斯 更 新 后 的 模型 参数 值 也 示 于 表 3 中 。MCMC 
分 析 时 选取 抽样 迭代 10 000 次 。 由 表 3 中 的 数据 可 
以 看 出 ,采用 MCMC 抽样 的 贝 叶 斯 更 新 模型 的 参数 
估计 值 的 标准 误差 均 比 表 2 中 未 更 新 模型 相应 的 参 
数 的 标准 误差 要 小 ,而 且 MCMC 分 析 和 迭代 轨迹 很 忆 
稳定 在 样本 的 均值 ,参数 的 自 相关 性 很 快 消失 ,这 证 
明 MCMC 分 析 过 程 是 稳定 收敛 的 ,抽样 结果 是 有 效 
的 。 另 外 通过 MCMC 标准 误差 ( MCSE) 的 数值 来 看 ， 
所 有 参数 估计 值 的 MCSE 数值 均 小 于 标注 误差 的 
5% ,这 也 说 明 通过 MCMC 抽样 算法 的 贝 叶 斯 更 新 广 
法 对 机 场 道 面 性 能 退化 过 程 的 威 布尔 分 布 生存 分 析 
模型 的 参数 估计 的 后 验 分 布 接近 于 其 真实 的 分 布 。 

一 下 文采 用 新 模型 的 参数 分 别 计算 出 机 场 道 面 性 
能 退化 过 程 的 5 个 状态 的 转移 概率 ,融合 新 旧 刚性 
遵 沁 数据 ,采用 1 年 的 检测 时 间 间隔 ,计算 出 机 场 道 
面 样本 均值 性 能 退化 过 程 的 转移 概率 矩阵 ,分 别 作 
名 态 1 ~ 状态 4 的 转移 概率 和 状态 持续 时 间 的 关 
系 图 ,如 图 1 ~ 图 2 所 示 。 为 了 方便 比较 ,将 新 \ 旧 模 
区 的 转移 概率 一 起 表示 在 图 中 ,其 中 * 己 "表示 新 模 
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(a) 状 态 1 评价 等 级 优 


图 1 机 场 道 面 性 能 从 状态 并 评价 等 级 ) 天 
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型 状态 i 转移 到 状态 j 的 概率 ,“P';” 表 示 旧 模型 状 
态 i 转移 到 状态 j 的 概率 。 

从 图 1 ~ 图 2 可 以 看 出 :新 模型 的 状态 转移 概率 
相 比 有 旧 模型 的 状态 转移 概率 的 变化 趋势 在 不 同 状 态 
等 级 表现 出 显著 的 不 同 ;状态 1 一 状态 2、 状 态 4 一 状 
态 5 的 新 旧 模 型 转移 概率 数值 变化 不 大 ,图 形 基本 
重合 ;状态 2 一 状态 3\ 状 态 3 一 状态 4 的 新 旧 模 型 转 
移 概率 数值 变化 较 大 ,图形 也 表现 出 明显 的 不 同 。 
这 些 和 新 加 入 的 机 场 道 面 PCI 数据 变化 特点 有 直接 
关系 ,新 加 入 的 机 场 道 面 性 能 衰变 过 程 对 于 机 场 道 
面 衰变 过 程 的 初始 状态 和 吸收 状态 影响 较 小 ,而 对 
于 机 场 道 面 性 能 衰变 过 程 的 中 间 状 态 则 有 和 较 大 影 
响 。 其 主要 原因 可 能 和 机 场 道 面 的 使 用 性 能 特点 有 
一 定 关系 ,新 加 入 的 机 场 道 面 样本 的 使 用 年 限 虽 然 
已 经 接近 使 用 寿命 ,但 是 由 于 机 场 运 维 部 门 对 道 面 
进行 日 常 养护 和 维修 ,及 时 解决 道 面 出 现 的 表 观 病 
害 和 结构 性 损伤 ,从 而 使 道 面 的 使 用 性 能 得 到 了 及 
时 恢复 ,促使 PCI 数值 也 逐渐 增 大 ,新 加 入 样本 的 性 
能 衰变 过 程 基本 均 处 于 中 间 状 态 ,进而 对 机 场 道 面 
性 能 衰变 过 程 的 中 间 状 态 产 生 了 较 大 影响 。 
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(b) 状 态 2 评价 等 级 良 
F 始 后 :年 的 转移 概率 (状态 1、 状态 2) 


Fig.1 The transition probability of airport pavement performance in t years after the beginning of state i 
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(a) 状 态 3 评价 等 级 中 


0 5 10 0 
持续 时 间 /年 


(pb) 状 态 4 评价 等 级 次 


图 2 机 场 道 面 性 能 从 状态 访 评 价 等 级 ) 开 始 后 1 年 的 转移 概率 (状态 3 状态 4) 


Fig.2 The transition probability of airport pavement performance in t years after the beginning of state i 


(evaluation grade) (state 3 ,state 4) 
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按照 第 1 节 所 述 的 方法 ,利用 贝 叶 斯 更 新 后 的 
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加 ; 同 理 ,对 于 数据 Data2S22 机 场 道 面 实 测 PCI 数 


参数 模型 和 半 马 尔 可 夫 随 机 过 程 , 求 出 机 场 道 面 任 
意 状态 等 级 在 任意 时 间 的 转移 概率 矩阵 P, (1) , 利 
用 转移 概率 矩阵 来 预测 机 场 道 面 在 任意 时 刻 ; 的 状 
态 或 等 级 。 为 了 对 比 新 旧 模 型 ,仍然 以 数据 集中 第 
8 组 数据 和 第 22 组 数据 作为 分 析 案例 ,得 到 该 道 面 
的 使 用 性 能 概率 预测 结果 。 新 模型 两 组 道 面 的 使 用 
性 能 预测 结果 如 图 4 所 示 。 

分 别 对 比 图 4(a) 和 图 3(a) .图 4(b) 和 图 3 
(b) ,可 以 看 出 ,新 模型 对 于 机 场 道 面 性 能 退化 过 程 
的 PCI 预测 曲线 更 符合 实测 的 PCI 变化 趋势 线 ,并 
且 新 模型 的 预测 区 间 基 本 都 包含 实测 的 PCI 数值 。 
从 数值 上 来 看 ,对 于 数据 DatalS08 机 场 道 面 实测 
PC 数值 衰变 到 50 的 时 间 为 24.7 年 , 旧 模 型 预测 的 
该 时 间 为 32 年 ,新 模型 预测 的 该 时 间 为 26.4 年 , 预 
洲 齐 误差 减 小 约 23% ,这 说 明 经 过 贝 时 斯 更 新 后 的 
模型 的 预测 精度 有 较 显著 提高 ,可 信和 度 有 较 显著 增 


©O 

《CO 

©O 

CD 

GN 

©O 

AN 

= 

>< 00 Ss 0 320” 23 30 35 0 4 “50 
© 持续 时 间 /年 
可 

ss (a) Datal S08 
© 


图 3 机 场 道 面 数据 样本 半 马 尔 可 夫 模 型 PCI 预测 图 


值 衰变 到 70 的 时 间 为 26. 8 年 , 旧 模 型 预测 的 该 时 
间 为 31.3 年 ,新 模型 预测 的 该 时 间 为 28.2 年 ,预测 
的 误差 减 小 约 12% 。 

对 于 预测 曲线 (PCI 值 状态 区 间 中 值 预 测 曲 线 ) 
在 早期 和 末期 均 偏 离 实测 PCI 数值 ,这 可 能 和 半 马 
尔 可 夫 过 程 的 无 后 效 性 有 一 定 关系 ,状态 持续 时 间 
随机 变量 的 概率 密度 分 布 不 仅 取决 于 当前 状态 ,而 
且 还 和 下 一 个 状态 相关 。 在 PCI 衰变 早期 ,跑道 的 
初始 PCI 数值 的 基本 上 都 是 100 ,所 以 实测 PCI 数值 
曲线 更 接近 于 PCI 预测 曲线 的 上 限 。 而 在 PCI 衰变 
末期 ,由 于 缺少 下 一 个 状态 的 PCI 数值 (因为 机 场 跑 
道 的 PCI 数值 一 般 小 于 50 时 , 就 需要 进行 维修 养 
护 ) ,再 加 上 机 场 跑道 日 常 维修 养护 的 作用 ,所 以 实 
测 PCI 数值 曲线 更 接近 于 状态 等 级 PCI 数值 区 间 的 
上 限 值 。 
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Fig.3 Predicted PCI relative to time instate( years) for airport pavement data samples with semi-Markov model 
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图 4 机场 道 面 数据 样本 半 马 尔 
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(b) Data2S22 


可 夫 模 型 PCI 预测 图 (新 模型 ) 


Fig.4 Predicted PCI relative to time in state( years) for airport pavement data 


samples with semi-Markov model( the new model) 


4 结论 


通过 对 两 种 机 场 道 面 性 能 检测 数据 的 融合 ， 
成 了 107 组 道 面 近 30 年 的 PCI 数据 来 估计 机 场 道 
面 性 能 退化 预测 模型 ,通过 分 析 , 结 论 如 下 。 

1) 机 场 飞行 交通 量 、 机 场 道 面 结构 面 层 混凝土 
厚度 以 及 机 场 道 面 投 入 使 用 年 限 对 退化 过 程 的 影响 
不 同 。 随 着 飞机 起 降 频次 的 增加 , 道 面 的 损坏 程度 
增 大 ,PCI 数值 减 小 ,机 场 道 面 维 持 当 前 状态 等 级 持 
续 时 间 的 危险 越 大 。 随 着 道 面 厚度 的 减 小 , 道 面 结 
构 承 受 答 裁 的 能 力 下 降 , 导 致 在 相同 飞行 交通 量 水 
平 下 道 面 的 损坏 程度 越 大 ,进一步 影响 PCI 值 的 减 
/人民 和 场 道 面 维持 当前 状态 等 级 持续 时 间 的 危险 


而 对 于 使 用 年 限 的 作用 ,除了 状态 等 级 为 良 的 
过程 模型 ,其 他 的 状态 等 级 均 没有 显著 的 统计 
学 意义 ,这 说 明 使 用 年 限 对 于 机 场 道 面 性 能 退化 过 
种 影响 作用 不 明晰 ,这 和 半 马 尔 可 夫 过 程 模型 的 无 
后 数 有 关 。 

2) 利 用 MCMC 进行 仿真 ,得 到 了 更 新 后 的 模型 


倩 评 结果 ,将 更 新 后 的 半 马 尔 可 夫 过 程 模型 的 预测 
结果 与 未 更 新 的 原始 模型 的 预测 结果 进行 比较 , 结 


ss 


明 采 用 MCMC 方法 的 贝 叶 斯 分 析 技术 是 一 种 
有 组 的 更 新 先 验 信息 的 反 分 析 方法 ,可 以 减 小 参数 
未 确定 性 , 较 显著 提高 预测 结果 的 置信 度 , 减 小 模 
更 卫 预 测 误差 。 

3 ) 通 过 随机 选取 数据 样本 进行 对 比分 析 , 结 果 
表明 ,更 新 后 模型 的 状态 转移 概率 相 比 原始 模型 的 
状态 转移 概率 的 变化 趋势 在 不 同 状态 等 级 表现 出 显 
著 不 同 ;新 加 入 的 机 场 道 面 性 能 训 变 数据 对 于 机 场 
道 面 误 变 过 程 的 初始 状态 和 结束 状态 影响 较 小 ,而 
对 于 机 场 道 面 性 能 衰变 过 程 的 中 间 状 态 则 有 较 大 
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